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The theory of nonlinear dipole plasma oscillations generated in a 
metal spheroidal nanoparticle by a laser-wave field has been devel- 
oped. Approximate (to within the cubic term) analytic expressions 
for the nanoparticle dipole moment have been obtained in the case 
where the laser field is oriented along the spheroid rotation axis. 



particle has one plasma resonance, whereas a spheroidal 
particle has two plasma resonances and an ellipsoidal 
particle has three ones. Therefore, the task aimed at 
developing the nonlinear theory of dipole plasma oscil- 
lations in asymmetric metal nanoparticles remains chal- 
lenging and interesting for today. 



1. Introduction 

When the center of masses of the electron subsystem in 
a metal nanoparticle is shifted with respect to the cen- 
ter of masses of the ion subsystem, there emerges an 
electrostatic force, which counteracts their spatial sepa- 
ration. This force may invoke dipole plasma oscillations 
in metal nanoparticles. In the first approximation, it is 
proportional to the relative displacement of the electron 
and ion centers of masses. If the displacement grows 
further, the electrostatic force starts to depend nonlin- 
early on this separation, which results in the appearance 
of nonlinear dipole plasma oscillations. These nonlinear 
plasma oscillations in metal nanoparticles were studied 
in works 0-01 • I n works 0,0], plasma oscillations were 
considered in the continual approximation, and the mi- 
croscopic approach was taken as a basis in work Q . In 
all cited works, the shape of a metal nanoparticle was 
assumed spherical. 

Our work is devoted to the development of the theory 
of nonlinear plasma oscillations in metal nanoparticles of 
ellipsoidal shape. It should be emphasized that the re- 
sults of the theory of plasma resonances in asymmetric 
metal nanoparticles cannot be reduced to small correc- 
tions to the results known for spherical particles, but 
have fundamental differences. In particular, already in 
the linear approximation, a spherically symmetric metal 



2. Formulation of the Problem 

We consider the problem of oscillations in metal 
nanoparticles in the continual approximation and take, 
as a basis, the hydrodynamic equations for the electron 
density n e (r, t) and the electron velocity v(r, t) similarly 
to work 0| : 



¥ + VM 



dv 
~dt 



(v V)v 



= 0, 

F 

m e 



(1) 



-eE L + e[V$ e + V$,] - — 



(2) 



In Eq. ©, F(r, t) is the total force that acts on the elec- 
tron liquid. It is composed of the action of the electric 
field generated by the laser wave, El, and the action 
of the gradients of electron, $ e , and ion, $i, potentials, 
and the pressure, p, gradient. In the dipole approxima- 
tion, the field Ei is considered spatially uniform within 
a nanoparticle. 

Let us introduce a vector that characterizes the posi- 
tion of the center of masses of the electron subsystem. 



u(t)=N~ 1 J d 3 rn e (r,t) . 



(3) 
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where iV e is the total number of electrons in the metal 
nanoparticle. With regard for Eq. (|TJ) , the equation of 
motion for the center of masses of the electron subsystem 
looks like 



„ (TU . , 



d 2 n e (r,t) 
dt 2 



(4) 



(5) 



/{ dv 
d 3 rn e (r,t) l— + (vV)v 

Rewriting it in the form 

m e u = JV- 1 (F), 

and using Eq. (|2J), we obtain 

(F) = —e E L (t)N e + 



+Jd 3 r {e [V$ e (r,t) + V$i(r,t)] n e (r,t) - Vp(r,t)}. 

(6) 

Above, we have briefly reproduced the approach to 
plasma dipole oscillations in metal nanoparticles used 
in work Q|. The authors of work [2], by making the 
required estimations, adopted an approximation, whose 
essence is the assumption that the electron subsystem of 
a nanoparticle shifts as a whole (without deformations) 
together with the center of masses of electrons. This 
allows us to put 

n e (r,t) = ni°)(|r-u(t)|), 



<S> e (r,t) = $(°)(|r-u(t)|), 



p(r,t)=po(|r-u(i)|). 



(7) 



We note that, at the thermodynamic equilibrium (i.e., 
in the absence of (£)), the equality 



(0) 







(8) 



is valid. In this case, Eqs. ([6]) and ([7]) yield 

(F) = -eE L (t)N e + J d 3 reV$f (r)n<°>(|r- u(t)|) = 



-eE £ (t) N e 



d_ 
du 



d 3 re^ 0) (r) nW(|r-u|). (9) 



(0)( 



We use this formula as the basic one. 

In contrast to the previous works where plasma 
nonlinear oscillations in spherically symmetric metal 
nanoparticles were considered, we analyze asymmetric 
nanoparticles. Let a metal nanoparticle have the shape 
of an ellipsoid of revolution (spheroid). Let the coordi- 
nate OZ axis be oriented along the spheroid symmetry 
axis. In addition, we suppose that the laser field Ei(t), 
as well as a shift of the center of masses of the electron 
subsystem induced by the field, is directed along the OZ 
axis. To emphasize this fact, we use the notation 



u = Z() 



(10) 



in what follows. It is worth noting that, if the particle's 
shape differs from the sphere, and the laser field that gen- 
erates plasma oscillations is not oriented strictly along 
the symmetry axis, the various plasma resonances are 
coupled with one another in the nonlinear approxima- 
tion, and the problem becomes incredibly complicated. 
Our calculations given below show that even the sim- 
plest model presented above and allowing for deviations 
from spherical symmetry is already capable to produce 
qualitatively new results, as compared with the spherical 
case. 

Up to now, except for formulas (0 and ((9]), the 
nanoparticle symmetry has not been specified. Similarly 
to work 0|, we adopt that the electron subsystem shifts 
as a whole (without deformations) together with its cen- 
ter of masses. Hence, we adopt that 

n e (r,t) = n o AoOR(0)-|r-zo|), (11) 

where no is the concentration of electrons, and Aq(x) is 
the step-like function 



A (x) 



1, x > 0, 
0, x<0. 



(12) 



The function R(6) defines the spheroid surface (to be 
more specific, let the spheroid be prolate), 



Ra 



R(9) = 



where e„ is the spheroid eccentricity, 



(13) 



R\ 
R 2 



- 1 



(14) 
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9 is the angle between the axis OZ and the R(0) on the 
spheroid surface, and and R± are the longitudinal 
(along the symmetry axis) and transverse, respectively, 
curvature radii. 

Supposing, in analogy with Eq. ([7|) , that the structures 
of functions $ e (r, t) and p (r, t) are similar to that of the 
function n e (r,t) given by formula (jllj) . we obtain 

(F) = -eE L (t) N e + n Q x 



x / AeV$f ) (r)A (R(8) - \r - z \) 



(15) 



instead of formula ©. Hence, in order to determine the 
force (F) that counteracts a displacement of the elec- 
tron subsystem in a spheroidal metal nanoparticle along 
the axis OZ, we have to determine the ionic electrostatic 
potential $|°' ) . This will be done in the next section. 

3. Electrostatic Potential of a Charged Spheroid 

The electrostatic potential generated by the ion core in 
a spheroidal nanoparticle looks like 



d 3 r 



r — r 



(16) 



Let the density of ion charges be uniformly distributed 
over the volume V , 



+2 Erri cos (^' ~ p n ( cos0 ') p n(^se). (20) 

z — ' (n + ml! 

The angles (tp',9') and (ip,8) in Eq. (2DD describe the 
spatial orientations of the vectors r' and r, respectively. 
When substituting Eq. (|2(jp in Eq. (|18l) and integrating 
the result obtained over tp' , the second term in Eq. (|2U)) 
is nulled. 

According to Eq. fig]), the quantity i?(6»') in Eq. §11 
satisfies the condition 



R± <R{0') < R\\ 



(21) 



Therefore, it is expedient to consider the integral over r' 
in three cases: 

a) r < Rj_, 

b) R±<r < R\\ , and 



c) Ru < r. 



(22) 



At r < in accordance with Eqs. (|T8|) - ([20|) . we ob- 
tain 

oo * 

$f } = 2tt Pi Pn (cos 0) / dO' sm9'P„ (cos6')x 



pi = e Ni Zi /V — const, 



(17) 



where Ni is the number of ions, and Zi is the charge mul- 
tiplicity. Taking the spheroid symmetry and the charge 
uniformity into account, Eq. (|16p takes the form 



2-7T 7T 



R(e') 



$' 0) = Pi J dip' J d& sin 9' J 



I „I2 



dr r 



r - r' 



(18) 



To carry out the integration in Eq. (fl8j) . it is expedient 
to make the expansion 



^ = f>(cos,){K0 
1 n=a K ~ KT 7 ) 



at r' < r, 
at r' > r, 



(19) 



where v is the angle between the vectors r and r', and 
apply the relation [4] 

P n (COS V) = P n (cost)') P n (cOS0) + 




(23) 



The further calculation of the integrals in formula (f2"B1 
has no difficulties. A more difficult situation arises in the 
case R± < r < R» , because, in accordance with Eq. (|2T|) . 
the variable r' changes in the same interval. Therefore, 
depending on the angle 0', the maximum value of r' can 
be both larger and smaller than r (see Fig. 1). It is 
expedient to introduce an angle 9\ , at which the ellipsoid 
and the sphere of radius r intersect; in other words, 



R(6') = R(9i) = r. 

In view of Eq. ([13]), Eq. (JUJI yields 

2\ >! V2 



cos 9\ = — i 1 

e„ 



r 



(24) 



(25) 
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Now, let us decompose the integral over 9' in Eq. (fTSj) 
as follows: 



d6' . 



d6'...+ / d0'...+ / d& 



(26) 



Figure 1 demonstrates that r' can be both larger and 
smaller than r in the intervals < 9' < 9\ and ■k — 9\ < 
9' < it. At the same time, r' is always smaller than r 
in the interval 9\ < 9' < 7r — 9\ . Substituting expansion 
(fT9|) in Eq. (p~8|) and dividing the integration interval in 
accordance with procedure (|26[) , we obtain 



f "/ 

$' 0) = P„ (cos 9)1 / dfl'sinf?' 

n=0 



P„(cos 9')x 




R(e') 



r 

L o 



dr'r 



Fig. 1 



7T/2! 



+ / d6" sin 9' [P n (cos 6") + P n (- cos 6>')] x 



+ / d& sin 6"P„ (cos 9')~ ( dr'r' 2 ( r - 



+ J d6' sin 9' P n (cos 6') x 

7T-01 



o 



(27) 



Making substitutions of the type 0' -r 7r — 0' in the last 
term of Eq. (|2"T|) , the whole expression (f2"T)) can be ex- 
pressed in the form 



f f 

$, (0) =27rp i £p„(cos0K / d6" 

n=0 I { 



smf x 



x [P„(cos6>') +P n (-cos0')] x 



R(e') 



dr 



R(S') 

X 1 - [ dr'r' 2 ( r - 
r J \r 
o 



, v n 



(28) 



One can see that all terms with odd powers of n disap- 
pear from sum (|28|) . Expression ([28]) describes ^j ' in 
the range R± < r < R» . 

At last, let us consider the case r > Pn , where r' < r. 
From Eqs. (TS]| and ([15]). we obtain 

oo * 

$(°) = 2tt,# ^ P„ (cos 6>) / d6»' sin0' x P„(cos6>')x 



1 



R(e') 



x- / dr r 
r J \ r 

o 



As is seen from expressions (|2Tfl) . 
write 

oo 

= 2 P n (cos0)*(n) - *(0)- 



+P 2 (cos 6>) * © + P 4 (cos 0) * Q 



(29) 



and (123) ■ we can 



(30) 
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r m (v) 



Fig. 2 

for the whole range of variation of the vector r. Integrat- 
ing in Eqs. (f2"3")l . {25}, and (12"§)) , we obtain the expressions 
for the coefficients *f?(n). In particular, at r < R±, we 
find 



r r 2 R 2 

$(0 = 2^ -- + -iln 
L 3 2e p 

using Eq. {23]). From Eq. 

r<R h 

(r 2 

*(0) = 2npA — (cos 0i - 



1 + e. 



; r < iij., (31) 



, we obtain that, at R± < 



.^1 



In 



1 + e p cos ^i 1 



1 — e p cos 0i 1 



2R\ 



COS 0i 



Ri<r< Ru 



At last, Eq. ([29ll implies that, at r > R\\, 



^(0) = ^PiR±R\\R\\/r, r>R {l 



(32) 



(33) 



Similar expressions for ^(2) are given in Appendix. 

As is seen from Eq. (|2"5)) . cos0i = at r = R±. There- 
fore, Eq. {32J coincides with Eq. ^ . At r = fly, we 
have cos 0i = 1, and expression (pj2"j) coincides with 

Passing from the ellipsoidal shape to the limiting case 



it is easy to see that \&(n) — >• for all n ^ 0. Con- 
cerning \&(0), Eqs. (|3"Tj) - (f3"3")l show that, at this limiting 
transition, 



*(0) ->■ = V i 



(3-^)/2R, at r < R 
i at r > i? 



(34) 



where V = is the sphere volume. Expression for 

in the form (f34"|) was used in [2j, while consider- 
ing nonlinear plasma oscillations in a spherical metal 
nanoparticle. 

4. Electrostatic Force 

Provided that the field E/,(t) is oriented along the axis 
0Z 7 and the particle shape has the adopted symmetry, 
the field (F) is also directed along the axis QZ, i.e., 

(F) = -eE L {t)N e + 



+n e 



Jd 3 rK V $[ 0) (r) A (i?(0) - | r - z 



(35) 



where Ko is a unit vector directed along the axis 0Z. 
The integration range over r in Eq. (|35| is defined by the 
condition that the argument in the step-like function Ao 
is larger than or equal to zero, i.e., 



R(6) -\r-z \ >0. 



(36) 



In the case where relation (|36[) is the equality, we obtain 
a root 



r = r m (v) = { i? 2 (0) - z 2 sin 2 v } 1/2 + z 



(37) 



As is seen from Fig. 2, the vector v m (v) corresponds 
to that point on the surface of a shifted spheroid, which 
is determined by the vector R(0) on the surface of the 
initial spheroid. At the shift, the vector R(0) moves in 
parallel to itself from point to point 0'. According to 
Fig. 2, we can write 



r m {v) com* = zq + R(9) cos0, 



R(0) sin0. 



(38) 



(39) 



Multiplying relations (|3"5)) and (|39p by smv and cos v, 
respectively, and subtracting the results, we obtain 



of the spherical shape (R± = Ru = R), i.e. at e p — )• 0, zq sinv = i?(0) sin(0 — v). 



(40) 
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This formula gives a relation between the angles 9 and 
v at a fixed zq. At zq — > 0, we obtain 9 — > v. Since the 
ratio zq/R{v) is small, we can write 



(41) 



so that Av can be determined from Eq. (|40|) by itera- 
tions: 



Aj/ 2 



A i/ 3 



R{v) 



R(u) 



sin v 1 



sin v 1 



R(v) 



Az/i 



Aj/ 2 



(42) 



To obtain the explicit dependence of r m (v) on the shift 
zq, let us expand Eq. (|37[) in a series (the function i?(#) 
should also be expanded with the use of relations (|41|) 
and 02])) 



r m (i/) =R(u)+z (l-e;)- 



cos ^ 



e 2 cos 2 v 



2 1 — e 2 sin 2 
Z ° 2i? ± (l-e 2 cos 2 i/) 3 / : 



- + 




(1 — e 2 cos 2 v)" 1 



cos ^ sm 



(43) 



For the spherical shape (e p —> 0), Eq. (j37|) yields 



r m (i/) = i?+z cosz/-^sin 2 ^-^sin 4 i/+... . (44) 

Comparing expressions (|4"3")l and (|4"4"1) . we see that, in the 
case of an asymmetric particle (e p ^ 0), an extra cubic 
nonlinearity absent in expression (|44l) emerges. 

To avoid a misunderstanding, we should emphasize 
that, in what follows, the term "asymmetric particle", 
will mean a particle, whose shape differs from the spher- 
ical one, rather than the absence of any symmetry ele- 
ments. 

Note that formula (|33|) with regard for Eqs. (|3U)) and 
(137[) can be written in the form 

(F) = -eE L (t)N e + 



.0) 



+27T7i e J dusmu 


Expressing r m (y) as 

r m (i/)=JZ(i/) + Ar(i/), 



dr r 2 



K V$f(r) 



(45) 



(46) 



and comparing it with Eq. (1431) . we see that only Ar(^) 
depends on the shift zq, and Ar -C i?(^). It would seem 
that the integral over r in Eq. (I45[) can be expanded in a 
power series in powers of Ar m (^) to obtain terms, both 
linear and nonlinear in zq. However, there is a subtle 
point. In order to obtain terms nonlinear in zq, we must 
differentiate the integrand in Eq. (|4"5)) . i.e. the function 



w(r) 



r 2 K V$f 



(47) 



However, as is seen, e.g., from the exact expression for 
the ion electrostatic potential in the spherical case (|3~4")) . 
the function ^"'(r) and its first derivative are continu- 
ous at the point r — R, i.e. across the particle surface. 
But already the second derivative of <&( s )(r) is discon- 
tinuous at this point. Therefore, the mentioned integral 
cannot be expanded in a Taylor series at the surface. 
At the same time, the ion electrostatic potential is de- 
scribed by smooth functions to the left and to the right 
from the surface. Therefore, we will do as follows. Let 
us divide the integration path over v in Eq. (|45|) into 
intervals, in which r m (y) lies only to the left or only to 
the right from the nanoparticle surface. Then, to the left 
and to the right from the surface, there exist the eligible 
reasons for function (|47|) to be expanded into a Taylor 
series. 

Let us explain the essence of our approach using, as 
an example, the spherical shape, for which the results 
are already known [l|,0|. Hence, the integral entering 
Eq. (|45|) can be rewritten as follows: 



dv sin? 



dr 1 



K V$ 



(0) 



tt/2 r m (v) n r m (v) 

= J dv sinw J drw(r)+ J dv sinw J drw(r). 

tt/2 

(48) 

As is seen from Eq. (|44l) . r m (v) > R for the first integral, 
and r m {v) < R for the second one (these inequalities 
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become somewhat violated at v 7r/2, but this does 
not make an appreciable contribution to the integral). 
In the case of the spherical shape in accordance with 
Eqs. ([34]) and ([3?]). we have 



(49) 



w(r) = —pi V cos u, at r > R, 



w(r) = —pt V cosv j at r < i?. 



(50) 



Now, let us substitute functions (|49|) and ([50]) into the 
first and the second integral, respectively, on the right- 
hand side of expression ([4"T]) . Then, let us take into con- 
sideration that r m {v) = R + Ar(y), with Ar(v) <C R, 
and make the expansion 



r m {v) 



dr r 



K V$f (r) 



_R+Ar(u) 

dr «;(r) 

o 



= / drw(r)+ w(r)Ar+- ( — ) (Ar) 2 
J 2\dr J R 



Taking into account Eqs. (|5Tj) and ([52"]) . all integrals in 
expression (|45[) can be calculated to the end, and we 
obtain (at zq > 0) 



(F) = - 



( 9 z 2 
e E L {t) N e - m e N e co 2 p i z - — 



32i? 3 



(54) 



where uj 2 = 47T 3 ( m n ° is the square of the plasma (dipole) 
frequency. 

Thus, we repeated the result of works Now, let 

us apply the same approach to the case of a spheroidal 
nanoparticle. 

5. Nonlinear Dipole Plasma Oscillations of 
Asymmetric Metal Nanoparticle 

In the case of a spheroidal nanoparticle with regard for 
Eq. ([30]), we have 



K V$ 



cosv 



(0) 



d<S> 



(0) 



sinv d<f> 



(0) 



cosu- 



dr 



dv 



fdtf(O) d*(2) 



\ dr 



ch 



P 2 (cosu)+ 



1 



d 2 ' 



1 



d 3 » 



24 V dr 3 



(51) 



To avoid the misunderstanding, we note once again that 
the derivatives in Eq. (|51| are not calculated at the sur- 
face (r = R), but at a point r, provided that r ap- 
proaches the surface from the left or from the right. In 
particular, in the given specific case, the matter concerns 
the expansion of function (|50p . provided that r — > R and 
r <R. 

From expressions (J43J) and ([44]) . we see that Ar is a 
power series in Zq, i.e. 



Ar = Ari + Ar 2 + Ar 3 + Ar 4 , 



(52) 



where An 



In particular, in the case of spherical 



symmetry, according to Eq. 



we have 



Ari = zqCosw, Ar 2 = 



~2R Sm W » 



Ar 3 =0, Ar 4 = -^smV 



-*(2) [l-P 2 (cosw)] 



(55) 



Below, we take into account explicitly only ^(0), al- 
though the required calculations were also carried out 
making allowance for the contribution of the function 
\£(2), the expression for which is presented in Appendix. 
Our estimates showed that the account of the contri- 
bution made by X P(2) does correct, to some extent, the 
coefficients at the powers of zo, but does not change our 
main conclusions. 

Hence, in accordance with Eq. ([4"3"]l , we have 

r m (0) = #|| +z > R ± ; r m {n) = R\\ - *o < R\\- (56) 

Taking these inequalities into account, let us split the 
integral over the angle v again as was done in Eq. (|48|l . 
We now substitute function (f33|) in the first integral on 
the right-hand side of Eq. (|4"5)k here, in accordance with 
Ea.lBl)! r m (0) > R\\. Then, according to Eqs. 07} and 
(|30| . we have 



(53) w(r) = cosv r 



2^(0) 



4?r 



dr 



Pi R± i?n cosv. 



(57) 
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In the second integral in Eq. (|48|) . for which r rn (n) < Ru , 
we also suppose that r rn (v) > R±_ and use function p2[) . 
In this case, we obtain 

w(r) = ^pi (l(r 2 -i?i) 3 / 2 -r 3 lcos«. (58) 
3 [e p J 

Note that the assumption r m (v) > R±_ also means that 



r(ic) — Ru — zq > R± or zq < R« 



Ra 



(59) 



Without assumption ((59]), the expressions for the ion 
electrostatic potential, as well as the integration limits 
(j4~4"]) , are transformed into the corresponding result for 
a spherical particle at e p — > 0. Condition (|59|) makes 
this passage to the limit impossible, because, if e p — > 0, 
inequality (|59|) becomes invalid at any small, but finite 
value of zq. In this case, for the passage to the limit 
e p — > to be eligible, one should engage function pip 
rather than function (j3"2"]) . Hence, the substitution of 
Eqs. ([57} and ([58]) in Eq. (08} gives 



sin v 



dr 



KnV$ 



(0) 



47T 9 
= __ J? 2 E||p . 



47T 



tt/2 



dvr m (v) sinw cosw + 



Pi J dvS ULVCO S V 
n/2 



(v) 

dr 



-{r 2 -Rlf 



(60) 



To integrate the second integral over r in Eq. (|60| . we 
take advantage, similarly to Eq. (f5Tjl . of the smallness 
of quantity Ar and expand this integral in a series in 
Ar. However, there exists a certain difference between 
cases (fBl~T) and (|60|) . In case (fBTj) . r TO (u) = R + Ar(v), 
and in case (|60|) . r m (u) = i?(w) + Ar(w) in accordance 
with Eq. (j43|) . To avoid excess complications, we ex- 
pand the integral Jq™^ dr {. . .} into a series in Ar at 
the point i?||, rather than at R(v). A reason for this 
approximation is that, first, the electrostatic potential 
is mainly governed by the distribution of charges near 
the ellipsoid vertex (pole), i.e. by the range of angles, 
in which R(v) w Rn, and, second, the function ^(0) and 
its derivative, as is seen from Eqs. (|32|) and ([33]), are 
continuous at the point r = Ru , similarly to what takes 
place in the spherical case, for which the exact solution 
is known. 



From Eq. (|60j) , confining the expansion to terms cubic 
in Ar, we obtain 



K V$ 



(0) 



r m (v) 

dv sin v J dr r 2 


tt/2 

= —V pi J dv Ar(v) sin v cos v- 
'0 



-V pi J dv sinu cosu 

TV /2 



Ar(v) 



(Ar(v)) 3 
2elR\ 



(61) 



where V = ^-R\R\\ is the volume of spheroid. 

Note that, owing to inequality (|S"9"]) , the following in- 
equality, as can be easily verified, is also valid: 



e p > Zq. 



(62) 



If inequality (|6"!?|) is obeyed, the term cubic in Ar is much 
smaller than the linear one, as it must be when expand- 
ing in a small parameter. Since the expression for Ar(v) 
itself is a series expansion in zo, we confine the consider- 
ation below to the terms, the order of which is not higher 
than Zq, i.e. we make the substitution 

Ar(v) » Ar 1 (v) + Ar 2 {v) + Ar 3 (v), 



(Ar(v)f w (A ri (v)f 



(63) 



into Eq. (l6"Tj) . The form of Ari(w)-terms for the 
spheroidal shape is clear from expression (|4"31 . 

Calculating the corresponding integrals in Eq. ([BT]) 
and substituting the obtained expression into Eq. (|45"]) . 
we obtain 



(F) = -e E L (t) N e - m e N e lo 2 z - m e N e u. 



S(e p ) 



2pL Rl 



No- 



where 



L\\u 2 pL 



2el 



In 



l + e ? 

1 - e. 



2e t 



u P Li 



(64) 



(65) 



^pL = *J is the plasma frequency, and L|| is the 

depolarization factor along the symmetry axis in the case 
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of prolate spheroid > R±) . In addition, we intro- 
duced a dimensionless parameter 8{e p ) in Eq. (|64p . which 
depends only on the eccentricity e p and looks like 



6(e p ) 



315 1 



13 

12if 



2e* 



5 - 6e 2 + el (\ + e, 
— - ' In 



8e p 




(66) 



Comparing the expressions obtained for the electrostatic 
force in the cases of spherical (formula (|54jl ) and ellip- 
soidal (formula (|6~4")0 nanoparticles, we see that, for the 
asymmetric particle, the quadratic nonlinearity changes 
to the cubic one. We recall once more that expression 
(|66p was obtained in the assumption < e„ < 1, and, 
therefore, the passage to the limit e p — > or e p — > 1 
cannot be justified. To get some understanding con- 
cerning the magnitude of parameter 8(e p ), we give the 
following values: 



In our case, when oscillations occur along the symmetry 
axis (u = Zo), the equation of motion ([5]) with regard for 
expression (|64|) reads 



(67) 



Zq + W|| z + io pL z — — . 



In Eq. (|67| . the dissipation was not taken into account. 
Therefore, solution has a singularity at to — » wii. 
The insertion of a standard term 2"/zq, which takes the 
oscillation attenuation into account, into the left-hand 
side of Eq. (|67j) gives rise to a disappearance of the sin- 
gularity from the solution. In particular, in the linear ap- 



proximation, instead of the solution 



eEp cos u t 



we have 



eE a 



cos(u t — <po) 



, where ipo is the phase. Similar 

modifications must also be made in the nonlinear terms. 
Note that Eq. J67|) with E L (t) in the form ([68]) corre- 
sponds to the so-called Duffing equation. The analysis 
of its solutions can be found, e.g., in work [5]. 

Now, having the explicit expression for the displace- 
ment of the center of masses of electrons, zq, in terms 
of the laser field Eq (see Eq. ([69])). we can write down 
the formula for the dipole moment of a spheroidal metal 
nanoparticle. If the laser field is polarized along the sym- 
metry axis of the spheroid, the dipole moment has the 
same orientation and equals 



d = Ve hq zq = di + d,2- 



(70) 



Here, d\ is the linear dipole component, which, in accor- 
dance with Eq. (|69|) . equals 



V 



dx = -r- 



47T uj - CJ 2 



Eq COSUjt. 



(71) 



Similarly, the cubic component of the dipole, (^3, can be 
written down in the form 

_ V5(e p ) (em) 2 ^ t 



The frequency u\\ corresponds to the frequency of a 
dipole plasmon, when the dipole oscillates along the sym- 
metry axis of the spheroid. If we put 



E L (t) = E cos cot 



(68) 



and assume that the nonlinearity is weak, Eq. (|67[) can 
be solved using the iteration method: 



zq 



eE coscot io^ L S(e p ) ( e E m e ) 3 



■ 00' 



4Ri 



3 cos co t 



cos 3lo t 
;' 2 - (3w) 2 



(69) 



3 cos to t 
w 2 — to 2 



cos Scot 



1 



(3w) 



At last, we would like to make the following remark. If 
the dimensionless displacement a = zq/ R\\ is introduced, 
Eq. (E7D looks like 



a + to 2 a + tOp L S(e p ) a 3 



eE(t) 
mRu 



(72) 



We see that the magnitude of nonlinearity is determined 
by the dimensionless parameter 8{e p ), the analytic form 
of which, as a function of the eccentricity e p , is given 
by formula (|66ll . In addition, the specific <5(e p )-values at 
e p = 1/2 and 1/5 were quoted above. This allows us to 
estimate the nonlinearity. 



1118 



ISSN 2071-0186. Ukr. J. Phys. 2011. Vol. 56, No. 10 



NONLINEAR PLASMA DIPOLE OSCILLATIONS 



It is also worth noting that the nonlinearity consid- 
ered above is induced by the electric component of the 
laser wave field. Under certain conditions (the particle 
size, the field frequency), the nonlinearity can be induced 
by the magnetic component. In particular, the effect of 
second harmonic generation in spherical metal particles 
under the influence of the magnetic component of the 
laser wave field was considered in work Q. 

6. Conclusions 

It has been shown that the laser field oriented along the 
symmetry axis of a spheroidal metal nanoparticle gener- 
ates a cubic nonlinearity, which is absent in the case of 
a spherical particle. Instead, the quadratic nonlinearity 
inherent to the case of spherical symmetry disappears. 
An approximate analytic expression for the dipole mo- 
ment of a spheroidal metal nanoparticle has been derived 
to within terms cubic in the field. 



47T 



= — PiR\< - J - 



15 



V r 



at r > Ru 
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APPENDIX 



*(2) = 2npir 2 
at r < R±. 



1 



(1 - cos 3 6»i) + - cos 2 6»i)cos6»i- 
4^(1-008 0!)- 



*(2) = 2rc Pi r 2 . 



■ In 



1 + ep cos 0\ 1 — e ? 
1 — e.p cos 6\ 1 + ej 



(3 c 



COS V\ 



npii R± < r < R\\ 



17. M. ToMHyK, fi.B. BymeuKo 
P e 3 K) m e 

y po6oii p03BHHyTO Teopiio HejilHifiHHX ^ihiojibhhx njia3MOBHx 
KOJiHBaHB y MeTajieBift HaHOMacTHHiji cc^epoiflajibHOi cJ^opMH, Hid 
reHepyiOTBCH nojieM jia3epHOi xBHjii. J\ji5l Buna^Ky jia3epHoro no- 
jih, opieHTOBaHoro b3,zj,ob:>k oci oSepTaHHH cc^epoifla, OTpmviaHO 
Ha6jiH^:eHi aHajiiTH^m Bnpa3H ^jih flHnojiBHoro MOMeHTy HaHO- 

MaCTHHKH (3 TOMHiCTK) ,0,0 Ky6iMHO! CKJia^OBOl) . 
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